Quenching of spectroscopic factors for proton removal in oxygen isotopes 
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We present microscopic coupled-cluster calculations of the spectroscopic factors for proton removal 
from the closed-shell oxygen isotopes 14 > 16 > 22 . 24 > 28 o with the chiral nucleon-nucleon interaction at 
next-to-next-to-next-to-leading order. We include coupling-to-continuum degrees of freedom by 
using a Hartree-Fock basis built from a Woods-Saxon single-particle basis. This basis treats bound 
and continuum states on an equal footing. We find a significant quenching of spectroscopic factors 
in the neutron-rich oxygen isotopes, pointing to enhanced many-body correlations induced by strong 
coupling to the scattering continuum above the neutron emission thresholds. 
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The concept of independent particle motion, and 
mean-field approaches based thereupon, has played and 
continues to play a fundamental role in studies of quan- 
tum mechanical many-particle systems. From a theoreti- 
cal standpoint, a single-particle (or quasiparticle) picture 
of states near the Fermi surface offers a good starting 
point for studies of systems with many interacting parti- 
cles. For example, the success of the nuclear shell model 
rests on the assumption that the wave functions used in 
nuclear structure studies can be approximated by Slater 
determinants built on various single-particle states. The 
nuclear shell model assumes thus that protons and neu- 
trons move as independent particles with given quantum 
numbers, subject to a mean field generated by all other 
nucleons. Deviations from such a picture have been in- 
terpreted as a possible measure of correlations. Indeed, 
correlations are expected to reveal important features of 
both the structure and the dynamics of a many-particle 
system beyond the mean-field picture. 

In a field like nuclear physics, where the average den- 
sity in nuclei is high and the interaction between nu- 
cleons is strong, correlations beyond the independent- 
particle motion are expected to play an important role 
in spectroscopic observables. In particular, experimen- 
tal programs in low-energy nuclear physics aim at ex- 
tracting information at the limits of stability of nuclear 
matter. Correlations which arise when moving towards 
either the proton or the neutron dripline should then pro- 
vide us with a better understanding of shell structure and 
single-particle properties of nuclei. So-called magic nu- 
clei are particularly important for a fundamental under- 
standing of single-particle states outside shell closures, 
with wide-ranging consequences spanning from our ba- 
sic understanding of nuclear structure to the synthesis 
of the elements P, Q ■ Unfortunately, the correlations in 
many-particle systems are very difficult to quantify ex- 
perimentally and to interpret theoretically. There are 
rather few observables from which clear information on 



correlations beyond an independent particle motion in a 
nuclear many-body environment can be extracted. 

A quantity which offers the possibility to study devia- 
tions from a single- particle picture, and thereby provide 
information on correlations, is the spectroscopic factor 
(SF). From a theoretical point of view they quantify what 
fraction of the full wave function can be interpreted as an 
independent single-particle or hole state on top of a corre- 
lated state, normally chosen to be a closed-shell nucleus. 
Although not being an experimental observable per se 0- 
5], the radial overlap functions, whose norm are the SFs, 
are required inputs to the theoretical models for nucleon 
capture, decay, transfer and knockout reactions. There is 
a wealth of experimental data and theoretical analysis of 
such reactions for stable nuclei[l|,@,0|. Data from (e, e'p) 
experiments on stable nuclei |l| indicate that proton ab- 
solute SFs are quenched considerably with respect to the 
independent particle model value, with short-range and 
tensor correlations assumed to be the main mechanism. 
Adding long-range correlations as well from excitations 
around the Fermi surface, one arrives at a quenching of 
30 — 40%, see for example Ref. Q. Nuclear physics offers 
therefore a unique possibility, via studies of quantities 
like SFs, to extract information about correlations be- 
yond mean-field in complicated, two-component, many- 
particle systems. 

Recent data on knockout reactions on nuclei with large 
neutron-proton asymmetries indicate that the nucleons of 
the deficient species, being more bound, show larger re- 
ductions of spectroscopic strength than the more weakly 
bound excess species 0, [1(3] • It is the aim of this work to 
understand which correlations are important when one 
moves towards more weakly bound systems. For this, we 
study the chain of oxygen isotopes and compute SFs for 
proton removal from 14 > 16 > 22 : 24 > 28 0. These isotopes span 
a large range of proton- neutron asymmetries, from 4/3 in 
14 to 2/5 in 28 0. Using ah initio coupled-cluster theory 
described below , we argue that the reduction in SFs 
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is due to many-body correlations arising from the cou- 
pling to the scattering continuum in neutron-rich oxygen 
isotopes. After these introductory remarks, we give a 
brief overview of our formalism, before presenting our 
results and conclusions. 

The spectroscopic factor -S^_ 1 (Zj) = jO^_ 1 (Zj;r)| , is 
the norm of the overlap function, 



Oi-M, r) = f (A - 1 1| a nlj || A) <f> nlj (r). 



(1) 



Here, O^_ 1 (0; r ) is the radial overlap function of the 
many-body wave functions for the two independent sys- 
tems with A and A — 1 particles, respectively. |^4) and 
\A — 1) can in general either be in their ground- or any 
excited state. In this work we consider only overlaps with 
\A) and \A — 1} in their ground states. The double bar 
denotes a reduced matrix element, and the integral-sum 
over n represents both the sum over the discrete spec- 
trum and an integral over the corresponding continuum 
part of the spectrum. The annihilation operator a n ij is 
a spherical tensor of rank j. The radial single-particle 
basis function is given by the term </> ra jj(r), where I and 
j denote the single-particle orbital and angular momen- 
tum, respectively, and n is the nodal quantum number. 
The isospin quantum number has been suppressed. We 
emphasize that the overlap function, and hence also its 
norm, is defined microscopically and independently of the 
single-particle basis. It is uniquely determined by the 
many-body wave functions \A) and \A— 1). From the 
definition of the overlap function in Eq. ([T]) it is clear 
that the SF is mainly a measure of how well nucleus A 
can be described by a single, uncorrelated nucleon at- 
tached to nucleus A — 1. Large deviations from unity 
indicate an increased role of many-body correlations be- 
yond a mean-field picture. For calculational details see 
Ref. [13 ■ 

We use the coupled-cluster (CC) ansatz [llj \ipo) = 
exp (T) |</>o) for the ground states of the closed-shell oxy- 
gen isotopes 14 > 16 > 22 > 24 > 28 0. The reference state, \cj>o), is an 
antisymmetric product state for all A nucleons. The clus- 
ter operator T introduces correlations as a linear combi- 
nation of particle-hole excitations T — T\ + T 2 + . . . + Ta, 
where T n represents an n-particle-n-hole excitation op- 
erator. For the CC singles and doubles approximation 
(CCSD) employed in this work, T is truncated at the 
level of double excitations, T = T\ + T 2 . 

Due to the non-hermiticity of the standard CC for- 
malism, we need to calculate both the left and the 
right eigenvectors. These are determined via the 
cquation-of-motion CC (EOM-CC) approach as \A) ss 
\R*{J A )) = exp(T)i?^(J A )|0 o ) and (A\ « (Lf(J A )\ = 
{(t>o\L^(J A )exp(-T). The operators R^(Ja) and 
L^(Ja) produce linear combinations of particle-hole ex- 
cited states when acting to the right and left, respec- 
tively. In the spherical form of the EOM-CC approach, 
the operators have well defined angular momentum by 



construction, as indicated by J A , which stands for the 
angular momentum considered. If the A-body system is 
in its ground state, the right EOM-CC wave function is 
identical to the CC ground state. 

Solutions for the A — 1-body systems are ob- 
tained with particle-removed equation-of-motion 
coupled-cluster method, where we use the CCSD 
ground state solution of the closed-shell nu- 
cleus A as the reference state in order to deter- 
mine the corresponding left and right eigenvectors 
L4-1) 
(A-l\ 



\Rj-\J A -i)) 



expCO^-V^-i)!^)) and 
= (^ol^-VA-OexpC-T). 



In actual calculations, the EOM-CC wave functions 
are obtained by determining the operators R^(J A ) and 
L^(J A ) as eigenvectors of the similarity-transformed 
Hamiltonian, H = exp (— T)i?exp (T). We refer the 



reader to Refs. 



for details about EOM-CC. 



Finally, we can approximate the SF in the spherical 
CC formalism as 



\Ja-i)\M\R${Ja)) 



x(R*-\j A ^)\\~a nl A\Lt(J A )) 



(2) 



where we have used the similarity-transformed spherical 
annihilation operator defined in Ref. 12]. The labels fi 
and v are included to distinguish between states in \A) 
and LA -1). 

The intrinsic A-nucleon Hamiltonian reads H — T — 
T cm + V, where T is the kinetic energy, T cm is the kinetic 
energy of the center-of-mass coordinate, and V is the 
two-body nucleon-nucleon (NN) interaction. We_employ 
here the N 3 LO model of Entem and Machleidt [15(. This 
interaction model is constructed with a cutoff of A = 500 
MeV. CC calculations starting from this Hamiltonian 
have been shown to generate solutions that are separa- 
ble into a Gaussian center-of-mass wave function and an 
intrinsic wave function, see for example Refs. 3, 16j |. 

We use a Hartree-Fock (HF) solution for the reference 
state, as detailed in for example Ref. ■ These HF 
solutions were built from the standard harmonic oscilla- 
tor (HO) basis combined with Woods-Saxon (WS) single- 
particle bound- and scattering states for selected partial 
waves. The role of the continuum is expected to be im- 
portant close to the dripline, as seen in Refs. 13, 17, lij |. 



For this purpose we use a spherical WS basis for the 
neutron S1/2, ^3/21 and d 5 / 2 partial waves. The single- 
particle bound and scattering states are obtained by 
diagonalizing a one-body Hamiltonian with a spherical 
Woods-Saxon potential defined on a discretized set of real 
momenta. We employ a total of 30 mesh points along the 
real momentum axis for each of the Si/ 2 , d 3 / 2 , and d 5 / 2 
neutron partial waves. For the harmonic oscillator basis 
we included all single-particle states spanned by 17 major 
oscillator shells. 

Figure [T] shows the calculated SFs for removing a pro- 
ton in the V\j2 and p 3 / 2 partial waves of 14 . 16 , 22 , 2 4,28q^ 
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FIG. 1: (Color online) Normalized spectroscopic factors 
for p\/2 and p 3 / 2 proton removal from the oxygen isotopes 

X4 lG 22 24 28 

' ' ' ' O. The continuum states included in the calcu- 
lation (HF-WS) lead to a dramatic quenching of the spec- 
troscopic factors as the neutron dripline is approached. For 
comparison, we show calculations of spectroscopic factors us- 
ing a HF basis built entirely from harmonic oscillator basis 
functions (HF-OSC). 



We compare our calculations of the SFs to calculations 
using an HF basis built entirely from harmonic oscilla- 
tor basis functions (HF-OSC, dashed lines). The results 
are obtained with an harmonic oscillator energy Hlj = 30 
MeV. Our calculations of the SFs depend weakly on the 
harmonic oscillator frequency, see for example Ref. 
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The P1/2 and p 3 / 2 proton orbitals are close to the Fermi 
level. In a traditional shell-model picture we would there- 
fore expect SFs close to unity for such states. However, 
we find a significant quenching of the SFs due to the 
coupling-to-continuum degrees of freedom. The calcu- 
lations done with a HF-OSC basis show no significant 
quenching, and illustrate clearly the limitation of the 
harmonic oscillator basis representation of weakly-bound, 
neutron-rich nuclei. This observation agrees also nicely 
with the analysis of Michel et al. [l9( . There, the authors 
demonstrate that the energy dependence of SFs due to an 
opening of a reaction channel can only be described prop- 
erly in shell-model calculations if correlations involving 
scattering states are treated properly. 

In our calculations the closed-shell oxygen isotopes 
14,16,22,24,28q are & rj bound with respect to neutron emis- 
sion (for this particular N 3 LO interaction with cutoff 
A = 500 MeV). In particular, we get 28 bound by 3.67 
MeV with respect to one-neutron emission. However, 
starting from a N 3 LO with a cutoff A = 600 MeV, we get 
28 unbound with respect to four-neutron emission and 
24 0, as seen in Ref. [20j. Clearly, three-nucleon forces 
are needed, to decide whether theory predicts a bound 
or unbound 28 0. We also computed SFs for the proton 
removal from 14 - 16 . 22 using the A = 600 MeV N 3 LO 



interaction model, and found similar results as for the 
A = 500 MeV N 3 LO interaction model. 

To further understand the role of correlations beyond 
mean-field we computed the SF for pi / 2 proton removal 
from 24 using the N 3 LO interaction evolved to a lower 
momentum cutoff A through similarity renormalization 
group methods [2l[. We used three different approxima- 
tions to \A) and \A — 1) and considered three different 
cutoffs A = 3.2, 3.4, 3.6fm _1 . First, in the crudest ap- 
proximation, using a mean- field HF solution for \A) and 
\A — 1), the SFs are by definition equal to unity. Sec- 
ondly, we used a HF solution for \A) while \A — 1) was 
approximated by one-hole and two-hole-one-particle ex- 
citations on the HF ground state \A), in this case we 
observed about 15 — 20% reduction in the SFs. Finally, 
our EOM-CC approach in Eq.f2|), gave a reduction of 
20 — 25% over the range of cutoffs considered. This 
clearly shows the importance of correlations beyond the 
mean-field. Changing the cutoff from A = 3.2£m _ to 
A = 3.6fm~\ the SF varied from 0.79 to 0.75, illustrat- 
ing the role of short range correlations on the quenching 
of the spectroscopic factors for proton removal in 24 O. 

The SFs of proton removal in the oxygen isotopes are 
determined by the squared norm of the overlap functions 
of Eq. JT]) . In order to probe the sensitivity of the tail of 
the overlap functions as we move towards 28 O, we com- 
pute the ratios of the absolute square of the radial overlap 
functions to the |( 15 N|a;j | 16 0)| 2 overlap function. These 
results are shown in Fig. [2] for the pi/ 2 proton state (the 
P3/2 proton state shows a very similar pattern). A no- 




FIG. 2: 



(Color online) Ratio of the radial overlap functions 
^ " . ;ii | 16 0),( 21 N|a !j | 22 0>,( 23 N|ay| 24 0>,and 
( 27 N|ajj | 28 0) for thepi/2 single-particle state. 



< 13 N|a ;j | 14 0},< 15 NK.| 



table reduction of these norms towards more neutron-rich 
nuclei is seen. The downward dip of the overlap ratios 
at larger radii comes from the fact the P1/2 proton or- 
bital become more and more bound as more neutrons 
are added to 16 0. For 14 the pi/ 2 proton is less bound 
with respect to 16 O, resulting in a bend upward. As the 
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neutron dripline is approached, the one-neutron emission 
thresholds for the oxygen isotopes and their neighbor- 
ing nitrogen isotopes are getting closer to the scattering 
threshold. Clearly, the tail of the wave functions will 
play a more important role as the outermost neutrons 
get closer to the scattering threshold. It is exactly this 
effect we observe in our calculations for the SFs for proton 
removal. Using a HF basis of purely harmonic oscillator 
wave functions, the density in the interior region of the 
nucleus is overestimated, while the density is shifted to- 
wards the tail when using a basis with correct asymptotic 
behavior. One should note that the nitrogen isotopes for 
a given neutron number are more loosely bound than 
their corresponding oxygen isotones, and this is the es- 
sential reason for the reduction. For 28 O and 27 N, no 
experimental values are available but if 28 O exists it will 
be very loosely bound and we may assume that 27 N is 
unbound. 
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FIG. 3: Plot of calculated SFs as functions of the difference 
between the calculated neutron and proton separation ener- 
gies. The results are for the single-particle states closest to 
the Fermi surface. For protons these are the px/2-states. 

Finally, we show in Fig. [3] the SFs of the proton and 
neutron states closest to the Fermi surface (for protons 
the pi/2-state), as a function of the difference between the 
computed proton and neutron separation energies. The 
results here agree excellently with similar interpretations 
made in Refs. 0, [T3j ■ One sees clearly an enhancement 
of correlations for the more strongly-bound, deficient nu- 
cleon species with increasing asymmetry. 

In conclusion, we have found a large quenching for the 
spectroscopic factors of the deeply-bound proton states 
near the Fermi surface in the neutron-rich oxygen iso- 
topes. This can be ascribed mainly to many-body corre- 
lations arising from a proper treatment of neutron scat- 
tering states. These results agree nicely with the math- 
ematical analysis performed by Michel et al (l9| . This 



result for the oxygen isotopes is similar to what has been 
inferred from neutron knockout reaction cross sections 
for deeply-bound neutron states near the Fermi surface in 
proton-rich sd-shell nuclei 0, [1(3] • Clearly, more work is 
needed to confirm the connection; experiments for proton 
knockout from oxygen should be undertaken and many- 
body calculations for proton-rich, heavy nuclei need to 
be carried out. 
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